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非 承 载 区 周 回 压力 槽 对 椭圆 轴承 
动力 学 性 能 的 影 啊 


杨 林 杰 , 陈 俊 豪 , 王 襄 元 ,张旭东 , 曾 志 刚 , 吴 鲁 纪 
(郑州 机 械 研 究 所 有 限 公 司 ,450052 郑州 ) 


摘 要 :为 探究 顶 圆 动 压 轴承 非 承 载 区 开 周 向 压力 油 楼 及 供 油 压力 对 轴承 -转子 系统 非 线 性 动力 学 
影响 ,本 研究 通过 改进 分 离 变 量 法 ,建立 了 适用 于 在 动 压 椭圆 轴承 非 承 载 瓦 开 楼 并 通 油 的 计算 模 
型 ,得 到 了 不 同 长 径 比 时 油膜 压力 随 偏心 率 的 变化 规律 ,分 析 了 开 周 向 油槽 宽度 与 供 油 压 力 对 轴 
承 - 转 子 动力 学 特性 的 影响 。 椭 圆 轴承 非 承载 区 只 开 构 不 供 压 力 油 对 系统 稳定 性 影响 较 小 , 增 大 开 
宰 寅 度 及 油槽 内 供 油 压 力 能 有 效 抑 制 转子 在 轴承 里 的 半 频 涡 动 成 分 ; 开 楷 宽度 为 承载 瓦 块 宽度 的 
0.3 倍 , 且 通 油 压 力 增 大 到 0.8 MPa 时 ,转子 振动 的 半 频 涡 动 幅 值 大 大 降低 ,转子 振动 幅 值 明显 降 
低 , 转 子 临 界 转 速 可 由 5 300 r/min 提高 至 8 400 r/min。 
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Effect of circumferential pressure oil groove in non-loading area 
on dynamic performance of elliptical bearing 


YANG Linjie,CHEN Junhao, WANG Jingyuan ,ZHANG Xudong,CAO Zhigang, WU Luji 


(Zhengzhou Research Institute of Mechanical Engineering Co. ,Ltd. ,450052 Zhengzhou, China) 


Abstract :Im order to explore the influence of the circumferential groove of the non-load bearing lobe of the 
dynamic pressure elliptical bearing and the oil supply pressure on the nonlinear dynamics of the bearing- 
rotor system , this paper uses the separation variable method to quickly solve the approximate analytical so- 
lution of the oil film pressure distribution on the elliptical bearing,and then obtains the law of the oil film 
pressure with the eccentricity under different aspect ratios. The influence of slot width and oil pressure on 
dynamic characteristics of bearing rotor is analyzed. In the non-loaded area of the elliptical bearing ,only 
slotting and not supplying pressure oil has little effect on the stability of the system. Increasing the slotted 
width and the oil pressure in the slot can effectively suppress the half frequency vibration of the rotor in the 
bearing. When the slot width is 0. 3 times of the lower pad width , and the oil pressure increases to 0. 8 
MPa ,the half-frequency eddy amplitude of rotor vibration is greatly reduced ,the amplitude of rotor vibra- 
tion is significantly reduced ,and the stable operation speed can be increased from 5 300 r/min to 8 400 Imin. 
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椭圆 滑动 轴承 作为 转子 的 重要 支撑 部 件 在 工业 
中 应 用 十 分 广泛 ,为 提高 其 性 能 ,通常 在 轴承 非 承 载 
瓦 内 表面 开 半 周 槽 ,这 样 可 以 提高 润滑 油 的 流通 性 
能 ,进而 起 到 降温 的 作用 ,相关 研究 开始 于 上 世纪 九 
十 年 代 。 朱 均等 针对 椭圆 轴承 在 上 瓦 开 周 向 本 
进行 了 大 量 理论 研究 ,确定 了 开 档 参数 对 油 温 变 化 

的 影响 规律 ,并 对 轴承 设计 进行 了 优化 ,大 大 降低 了 
轴承 的 能 耗 。 但 开 槽 也 存在 缺点 , 王 龙 刚 “ 通 过 研 
究 油 槽 的 径 向 滑动 轴承 油膜 压力 分 布 情况 ,发现 当 
开 覃 深度 达到 一 定数 值 后 ,可 将 覃 的 开设 位 置 处 压 
力 视 为 零 ,轴承 动 压 对 其 影响 可 忽略 。 

过 吴 兆 景 等 5 通过 求解 开 槽 轴瓦 的 轴承 动 压 润 滑 
放生 组 ,和 得 出 开 权 后 轴瓦 的 承载 能 力 下 降 等 结论 。 
注解 决 这 个 问题 , 李 建 鸿 等 四 通 过 合理 地 布置 轴 向 
矩形 止 槽 并 优化 凹 本 参 数 的 方法 ， 大 幅度 提升 了 轴 
承 的 承载 和 E 力 ,降低 了 轴承 摩擦 能 耗 。 丁 伟 '" 通过 
研究 发 现 改变 单 棉 泣 滑 楼 在 滑动 轴承 上 位 置 , 不 
仅 器 响 承 载 力 大 小 ,还 会 影响 油膜 的 最 高 压力 值 与 
给 使 椭圆 轴承 开 覃 方式 得 到 了 推广 。 

局 随 着 应 应 用 中 转子 转速 的 提高 ,在 支撑 轴承 内 会 
/ii 枢 训 动 和 油膜 振 区 造成 转子 失 稳 ,这 些 现象 
运用 传统 的 刚度 阻尼 线性 理论 上 已 不 能 解析 ,需要 
利用 非 线 1 性 理论 来 预测 。 近 年 来 , 随 着 技术 的 发 展 ， 
关 津 非 线 性 问题 计算 方法 的 研究 越 来 越 多 。 目 前 常 
用 前 为 有 限 差分 法 与 有 限 元 法 "|, 但 这 2 种 方法 
在 育 足 精度 的 要 求 下 ,计算 速度 非常 慢 ,无 法 及 时 快 
速 地 分 析 转 子 在 轴承 中 的 运动 规律 。 

为 能 快速 得 到 轴承 油膜 力 , 短 轴承 与 长 轴承 理 
论 被 提出 ,但 两 者 所 得 出 的 油膜 力 与 实际 油膜 
力 相 差 黄 大 ,不 能 满足 目前 精度 需求 。 郑 铁 生 等 
通过 变 分 途径 将 雷诺 边界 条 件 转 化 为 一 个 参数 优化 
来 模拟 油膜 破裂 ,从 而 提出 轴承 非 线 性 油膜 力 近 似 
解 。 王 永亮 等 "在 动态 5 油膜 胃 界 上 ,提出 以 分 

离 谈 量 法 来 求解 雷诺 方程 ,得 到 圆 瓦 与 椭圆 轴承 非 
线性 油膜 力 近似 解析 解 。 李 锋 等 … 在 分 离 变量 法 
的 基础 上 ,将 解析 解 分 解 为 特 解 和 通 解 之 和 ,根据 轴 
向 压 差 的 边界 条 件 得 到 圆 瓦 轴承 的 近似 解析 解 。 这 
些 算法 上 的 突破 为 椭圆 轴承 稳定 性 分 析 提 供 了 有 效 
手段 ， 但 都 只 适用 于 轴 向 两 侧 无 压 差 且 压力 为 零 的 

径 向 滑动 轴承 模型 , 即 只 适用 于 圆 瓦 轴承 与 浮动 环 
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密封 的 情况 。 耿 欢 '” 探究 微 织 构 对 轴承 润滑 性 能 
的 影响 ,发 现在 动 压 极 值 处 布置 微 织 构 可 提高 轴承 

能 力 。 张 乾 龙 等 "通过 计算 确定 最 优 的 油膜 参数 ， 
再 反 推 椭圆 轴瓦 具体 的 几何 尺寸 ,保证 轴承 提供 足 
够 的 刚度 与 阻尼 。 在 实际 生产 中 , 圆 瓦 轴承 应 用 较 
少 , 多 数 为 椭圆 轴承 与 多 油 棉 轴承 等 ,因此 对 实际 应 
用 意义 不 大 。 与 此 同时 ,转子 动力 学 也 在 不 断 发 展 ， 
肖 志 富 等 ”发 现 转子 系统 的 偏心 影响 着 刚性 支承 
转子 系统 的 临界 转速 分 岔 模式 。 本 研究 在 上 述 研 
究 的 基础 上 ,利用 分 离 变 量 法 ,提出 将 方程 的 解 可 拆 
分 为 特 解 和 通 解 两 部 分 ,其 中 特 解 满足 轴 疝 压 差 分 
布 , 通 解 满足 方程 分 离 变 量 近似 解 ,二 者 之 和 为 实际 
油膜 压力 解 。 对 椭圆 轴承 非 承 载 瓦 开 压 力 槽 后 的 情 
况 进行 模拟 研究 ,得 到 椭圆 轴 承 在 轴 向 压 差 和 主动 
润滑 共同 作用 下 的 压力 分 布 ,拓展 了 该 方法 的 应 用 


范围 。 本 研究 结论 可 为 椭圆 轴承 在 高 转速 系统 中 的 
应 用 奠定 理论 基础 。 
1 椭圆 轴承 油膜 力 计 算 模 型 
1.1 Reynolds 控制 方程 
本 研究 的 动 压 椭 圆 轴承 非 承载 瓦 开 权 结构 示意 


图 如 图 1 所 示 , 非 承载 瓦 上 内 侧面 开 全 周 向 槽 ,通过 

槽 上 所 开 孔 通 压 力 油 。 根 据 研 究 需 要 做 以 下 假设 : 

润滑 油 为 不 可 压缩 流体 ,密度 不 随 压力 与 温度 变化 ， 

忽略 流体 惯性 的 影响 ,在 油膜 厚度 方向 压力 大 小 不 
发 生变 化 ,得 到 的 无 量 纲 Reynolds 方程 为 

9/moPY 9/moPY_o8 ,oH 

a | eel | 00 . OT Cn 

无 量 纲 化 参数 设置 为 P=p/[6n0 (R/C) ]; 

t=2/R;e = EE/C;X = WC = Y= 三 


— esing; 
ecosp; Y=%/( CO) ;Y= 7/(CO);H=h/C;r = 102; 
NA,=L/D; A=L,/D, 
其 中 :h 为 油膜 厚度 ;p 为 油膜 压力 ;pn 为 动力 和 儿 度 ; 
en z 为 轴承 长 度 方 向 坐标 ;0 为 转 
速 ;L, 为 上 瓦 瓦 块 半 宽 度 ;E 为 转子 偏心 距 ;R 为 轴 
i 为 转子 半径 ;C =R-r 为 椭圆 轴承 半 
径 间 际 。 
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1 非 承 载 瓦 开 槽 结构 示意 图 


Fig.1 Schematic diagram of elliptical bearing 


structure parameters 
1.2 ”油膜 压 力 计算 推导 


油膜 厚度 计算 方程 为 
H=1 +ecos(0 -09) (2) 
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2 (了 +2X)sing + (X -2Y) cos0 
(3) 
由 (0) =0 可 得 动态 油膜 边界 为 
0, .a Tsign( 全) + 
Y +2 y+2X 


Dsign(2Y —X) 半生 (4) 


"对 于 椭圆 瓦 来 说 ,由 在 全 局 坐标 (轴承 中 心 坐 
,下 的 无 量 纲 量 转化 成 以 轴瓦 曲率 中 心 0| 与 
坐标 下 的 无 量 纲 量 , 则 扰动 位 移 与 扰动 速度 分 
别 蕊 变 为 
所 上 瓦 
Xx/C,Y, =y/C+m,X, =%/(CO),Y, =y/(C0) 

下 瓦 
X=x/C,Y, =y/C-m, X=%/(C0),Y,=7y/(C0) 
式 中 ,m=6/C 为 轴承 椭圆 度 。 

由 于 上 、 下 瓦 块 的 曲率 中 心 不 在 同一 中 心 点 上 ， 
需 根据 瓦 块 的 自然 边界 条 件 与 动态 油膜 边界 将 上 瓦 
在 周 向 圆周 方向 划分 成 5 个 区 域 1-V ,如 图 2(a) 所 
示 , 求 解 以 瓦 块 曲率 中 心 为 坐标 原点 的 油膜 正 压 起 
始 角 。 

定义 0 为 油膜 压力 通 解 的 正 压力 起 始 角 ,0, 油 
膜 压力 通 解 正 压 力 终止 角 。9; 为 油膜 压力 通 解 负 压 
起 始 角 ,0, 油膜 压力 通 解 负 压 终止 角 和 。 根 据 坐 标的 
转化 ,代入 求解 油膜 正 压 起 始 角 度 公 式 (4) ,确定 各 
轴瓦 上 动 压 油膜 区 间 。 
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上 瓦 圆周 模 


(a) 上 瓦 结构 示意 图 


(0) 下 瓦 结构 示意 图 
2 椭圆 轴承 结构 参数 示意 图 


Fig.2 Schematic diagram of elliptical 

bearing structure parameters 
1) 上 瓦 边界 条 件 确 定 
当 0. 在 工区 域 时 ,9 =0,,0, =7/2 一 y; 
当 9, 在 工区 域 时 ,上 瓦 无 正 压力 ; 
当 0. 在 亚 区 域 时 ,9 =3m/2 +y,0, =0,+T; 
当 0,. 在 区 域 时 ,9 =3mX2 +y,9 =5m/2 -yi 
当 0, 在 V 区域 时 ,0, =0.,0, =Sm/X2 -y。 
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上 瓦 负 压 边界 条 件 为 

当 0 在 I 区域 时 ,0; =3mX2 +y,0, =2m +0,; 

当 0. 在 工区 域 时 ,0; =3mX2 +y ,0 =5m/2 -Yi; 

当 0. 在 亚 区 域 时 ,0; =0, + 三 ,0 =5m/2 -yi 

当 0, 在 区域 时 ,上 瓦 无 负 压力 ; 

当 0. 在 V 区 域 时 ,9 =Sm/2 +y,0 =0.+2T。 

2) 下 瓦 的 边界 条 件 确定 

下 瓦 正 压 边界 条 件 为 

当 0, 在 I 区 域 时 ,0, =w/2+y,0, =0. +TT; 

当 0 在 卫 区 域 时 , 9, =7/2 +y,0, =37/2 -yy; 

当 9, 在 五 区 域 时 ,0, =0,,0, =37/2 -yi 

当 9, 在 人 V 区 域 时 ,下 瓦 无 正 压力 ; 

当 9, 在 V 区 域 时 ,0, =Sm/2 +y,6 =0,+T。 
其 中 : 工 为 上 瓦 油 权 宽 度 ; 为 轴瓦 包 角 ;y = 
(180° -nn)/2;9 为 转子 相对 于 轴承 中 心 的 状态 角 。 
再 当下 瓦 未 开 周 向 槽 ,无 需 定义 负 压 边界 条 件 。 

SS 性 无 量 纲 雷 诺 方程 进一步 整理 可 得 


3 这 60.6) 3esin(0 -9) .oP(0,€) ,FP(0,é) _h(O) 
Cop’ H 90 oe F 
GN (5) 
去 在 枯 圆 轴承 上 瓦 中 心 位 置 上 开 圆周 槽 ,并 在 村 
内 给 润滑 油 一 定 的 压力 p, ,轴承 两 端 为 泄 油 端 ,压力 

浙 P =0。 无 量 纲 化 得 到 P, =pi/[6n0 (R/C)”]。 
“轴瓦 开 槽 深度 为 5 mm 与 最 大 油膜 厚度 相 比 多 
出 至 个 量 级 ,因此 覃 内 轴承 油膜 压力 受 转子 动 压 的 

影 澳 很 小 ,几乎 无 变化 ,因此 假设 槽 内 的 润滑 油 压力 

各 媒 都 相等 , 且 不 随 权 内 周 向 位 置 的 变化 而 改变 。 
同时 轴瓦 的 交界 处 为 轴瓦 自然 边界 条 件 , 油 膜 压力 

为 (0) 确定 了 方程 (5) 的 边界 条 件 。 

根据 该 方程 为 二 阶 线性 非 齐 次 椭圆 方程 ,无 法 

得 到 解析 解 ,因此 假设 方程 解 可 近似 拆 分 成 轴承 主 

动 润滑 解 与 轴 向 压 差 解 两 部 分 , 即 该 假设 下 的 特 解 

与 通 解 相 加 ,其 中 通 解 满足 一 般 的 轴瓦 无 轴 向 压 差 

润滑 , 特 解 满足 轴承 轴 向 压力 分 布 。 

P,(0,é) =P,(0,é) +P,(0,é) (6) 
特 解 为 


P,(0,€) =A(é) (7) 
通 解 为 
P 05) =M(é) .+ N(0) (8) 
代入 式 (5) 可 得 
9A(é) oN(0) .3esin(0 -9). 
ge SP M(é) 
oN(0) oOM(éE) CO) 
a +N(0) 起 评 =0 (9) 


则 特 解 满足 方程 可 得 

4(C6) _ 

3 =0 (10) 
可 以 看 出 轴 向 压 差 只 影响 轴承 轴 向 压力 分 布 。 
通 解 方 程 满 足下 列 方程 。 


oN(0) ,3gsin(0 ~-p) .oN(0) 
Mo) pe) a 
OM(E) Ek(O) 
N(0) ge 0 (11) 


由 于 拆 分 成 两 个 方程 ,方程 的 边界 条 件 也 随 之 
变化 ,其 中 通 解 满足 边界 条 件 为 式 (12) , 特 解 满足 
边界 条 件 为 式 (13)。 

P.(0, -A)=P, 

P.(0,A) =P, 

A( -A)=P, 

A(A) = 已 
P,.(0, -A) =N(0) 
P,(0,A) =N(0) 
M( +A) =0 (13) 
N(0)=N(0,)=0 
N(0;) =N(0,) =0 

由 特 解 可 求 得 轴 向 压力 方程 解 为 
4(E) =ciE+c， (14) 

根据 边界 条 件 式 (12) 可 计算 得 出 c,、c,。 确 定 

方程 解 为 


(12) 


Pi-P:, Pi+P; i 
根据 边界 条 件 式 (13) 求 解 轴承 主动 润滑 油膜 
正 压 区 。 文 献 [15 ] 已 明确 推导 出 求解 过 程 ,这 里 不 
再 歼 述 。 
通 解 与 轴 问 压 差 解 和 为 
P, a P, +P, 
P,(0,é) 二 2 入 €+ 2 过 
(es +e ~“-e*—-e *) i 


(e™ +e ™) 


(Y+2X) 2 (X-2Y) , 


(—Y-2X)W,(0) + (X-2Y)W,(0) - 
i i (16) 
(Y+2X)W, - (X-2Y)W, 


[ Wi(01) + 


(—Y-2X)W,(0) + (X-2Y)W,(0)] 
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式 中 :gs (01,0,),ée( -A,A); W (0) = | wae: 


W,(0) = | S00d0;W.(0) = {SadosW, = W,( 0,) - 


W( 0,) ;] =1,2,3。 
同 理 , 可 求 得 轴承 润滑 的 负 压 解 与 轴 向 压 差 解 
和 为 


P,-P, Pi+tP', 
BE 
(es+e ~“—-e" —e *) 

(e+e ™) 


人 机 


FY -2X)W,(0) + (X-2Y)W,(0) - 
> | (17) 
€ 2 LA 


a + (X-2Y) W(0,)] 
戒 昌 :ee (0,,0,) ,Ee ( -和 A, 和 A)。 其 解 与 上 述 的 通 
解 三 致 , 轴 向 区 间 则 为 &e ( -和 ,1) ,这 里 不 再 重复 
八 韦 。 即 得 到 在 非 承载 区 轴瓦 圆周 槽 内 通 压 力 油 后 
钙 避 完整 的 油膜 分 布 情况 。 

1.3> 油膜 力 的 计算 

> 对 油膜 压力 积分 可 得 到 油膜 力 。 


| F, 6 m1/2 sj 
人 -六 fs pcos Roos 


一 cosO 


_ 6u{ 2R” 2 A [ sing 
3 | | (0,é) doe 
_ GUOR- 站 (18) 
C fo 
在 积分 时 需 考 虑 油膜 不 能 承受 负 压 , 即 当 己 <0 
时 ,认为 油膜 破裂 ,破裂 处 的 油膜 压力 为 0。 
F=f Ef (19) 
当 椭 圆 轴 承 的 上 瓦 开 槽 并 有 供 油 压力 时 ,椭圆 
轴承 总 的 油膜 力 为 
A x F ux F lx F gx 
= + + (20) 
(a 人 | | pa 
式 中 :为 槽 内 供 油 压力 所 产生 的 力 ;f 为 上 瓦 轴 
承 油膜 力 ;为 下 瓦 轴承 油膜 力 。 
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2 计算 模型 对 比分 析 


2.1 验证 油膜 压力 分 布 


本 研究 的 近似 解 可 由 式 (16)、 式 (17) 计 算得 
到 。 计 算 使 用 参数 值 如 下 :X=0.5;Y=0. 5; 槽 的 宽 
度 工 =0.3L;A,=0.9;D=100 mm;C =0.2 mm; 轴 
瓦 包 角 ng 为 150°;m=0.2;p, =1 MPa;p, =0 MPa; 
从 =0.029 8 Pa. s;w =100 rad/s。 

为 验证 所 求 得 的 近似 解 ,利用 经 典 有 限 差分 法 
再 次 求解 方程 ,具体 计算 参数 为 :上 、 下 两 个 轴瓦 在 
圆周 方向 的 网 格 数 为 150 个 , 轴 向 网 格 数 为 100 个 。 
采用 等 距 差分 格式 将 雷诺 方程 进行 离散 化 , 单 瓦 轴 
承 的 网 格 总 数 为 15 000 个 ,油膜 压力 收敛 准则 为 

二 IP,,-P,,| -7 
名 疙 | 于 [ 一 和 10 (21) 

式 中 :P, 为 某 一 近代 过 程 中 节点 (i,j) 的 压力 值 ;P;， 
为 上 一 次 迭代 的 压力 值 。 图 3 为 近似 解析 解 与 有 限 
差分 法 求 得 的 轴承 上 瓦 与 下 瓦 无 量 纲 油膜 压力 分 布 
图 。 选 择 计算 的 转子 中 心 位 置 为 (0.5,0.5) ,状态 
角 为 135° ,最 小 油膜 厚度 处 位 于 上 瓦 , 因 此 上 瓦 的 最 
大 油膜 压力 分 布 要 大 于 下 瓦 的 。 

定性 上 看 ,油膜 压力 分 布 一 致 ,定量 上 其 最 大 无 
量 纲 油膜 压力 分 别 为 1.597 与 1.591, 最 大 偏差 仅 为 
0.38% 。 本 研究 所 提 近 似 解析 解 模型 计算 用 时 总 时 
长 约 为 0. 08 s, 有限 差分 法 计算 用 时 总 时 长 约 为 
48 s, 近 似 解析 解 的 计算 速率 是 差分 法 的 600 倍 。 可 
以 说 明 运用 本 研究 改进 的 计算 方法 可 以 在 保证 计算 
精度 的 前 提 下 提高 计算 效率 。 


2.2 验证 油膜 力 计算 结果 


为 了 进一步 验证 近似 解 的 计算 方法 ,再 分 别 运 
用 两 种 计算 椭圆 轴承 油膜 力 并 对 比 计算 结果 ,分 析 
在 不 同 长 径 比 的 情况 下 ,两 者 求 得 的 油膜 力 随 偏 心 
率 的 变化 规律 。 固 定 转 子 中 心 与 轴承 中 心 的 状态 角 
为 1353。, 预 负 荷 系数 为 0.2 , 覃 内 供 油 压 力 为 1 MPa， 
结果 如 图 4 所 示 ,红色 线条 为 有 限 元 法 计算 结果 , 蓝 
色 线 条 代表 本 研究 近似 解 。 可 见 ,在 较 大 长 径 比 
(0.1 ~2) 范 围 内 ,二 者 的 计算 结果 基本 一 致 。 

由 于 预 负荷 系数 的 存在 ,转子 在 轴承 内 的 最 大 
偏心 率 不 会 超过 1, 因此 计算 时 偏心 率 取 0 ~ 0.83。 
随 着 偏心 率 的 增 大 ,转子 中 心 与 轴瓦 中 心 的 状态 角 
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夹 角 大 小 随 之 变化 ,油膜 正 压 起 始 角 度 也 随 之 变化 。 
在 横 内 供 油 压力 与 轴承 主动 润滑 油膜 力 的 共同 作用 


45 65 
重 纲 4 53 351 > 
一 全 村 周 向 角度 /() 


(a) 近似 解析 解 求 得 上 瓦 油膜 压力 分 布 图 


sl 2 BS 
-全 蒜 周 向 角度 /( 〇 ) 


(c) 有 限 差分 法 求 得 上 瓦 油膜 压力 分 布 区 


202212.00120v1 


综 上 所 述 ,通过 本 研究 所 推导 的 计算 模型 与 有 
限 比 分 模型 计算 结果 的 对 比 可 知 ,本 研究 所 推导 的 
计 竹 模型 更 适用 于 快速 求解 非 线性 轴承 -转子 动力 


10° 


aXiv 


-+ 4-0.1 有 限 差分 法 
一 9 一 4=0.1 本 研究 近似 解 多 
-x j=0.5 有 限 差分 法 ， 
六 -0 5 本 研究 近似 解 三 


区 
2 


-一 全 1 有 限 差分 法 
-…A… A=1 本 研究 近似 解 
一 一 4=2 有 限 差分 法 

0 4=2 本 研究 近似 解 


无 量 纲 油膜 力 
三 


本 


102 
Ge GL V2 D3 4 VS U6 RT VB 9 
偏心 率 


4 不 同 长 径 比 下 两 种 模型 油膜 力 随 偏心 率 变化 规 得 


Fig.4 Variation of oil film force with eccentricity of the 


I 
I 


two models under different aspect ratios 
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下 ,油膜 力 在 一 定 范围 内 并 不 随 偏心 率 的 增 大 而 增 
大 。 但 随 着 长 径 比 的 增 大 ,无 量 纲 油膜 力 随 之 增 大 。 


> 
一 
un 


无 量 纲 油膜 压力 
一 
三 


党 
SS 
吕 


将 125 145 165 185 205 225 245 
周 向 角度 /(9) 


(b) 近似 解析 解 求 得 下 瓦 油膜 压力 分 布 攻 


5 165 185 205 225 2 
全 周 向 角度 /() 
(d) 有 限 差分 法 求 得 下 瓦 油膜 压 力 分 布 图 


图 3 近似 解 与 差分 法 求 得 压力 分 布 图 


Fig.3 Approximate solution and difference method to obtain pressure distribution diagram 


表 1 不 同 偏心 率 下 有 限 差 分 法 与 近似 解 
法 求 得 的 油膜 力 (A =1) 
Tab.1 Oil film force obtained by finite difference method and 


approximate solution method under different eccentricity 


偏心 率 0.1 0.3 0.5 0.7 0.8 


有 限 差分 法 0.62 0.50 0.41 0.91 3.532 
近似 解法 0.62 全 .S1 0.44 0. 89 3.21 
从 表 1 中 可 以 看 出 在 长 径 比 为 1 时 ,在 大 偏心 
率 为 0.8 时 两 种 方法 计算 得 到 的 无 量 纲 油膜 力 最 大 
误差 为 8. 8% 。 进 一 步 验证 了 近似 解析 解 的 准 
确 性 。 


3 动力 学 分 析 


本 部 分 探究 在 椭圆 滑动 轴承 非 承 载 瓦 上 开 不 同 
宽度 的 周 向 槽 与 不 同 供 油 压力 对 轴承 -转子 系统 动 
力学 性 能 的 影响 。 
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影 1 f 度 的 增 

3.1 动力 学 方程 向 。 由 图 可 知 ,转速 一 定时 , 随 着 开 槽 宽度 的 


假设 转子 为 刚性 ,质量 为 2M,e 为 转子 质心 偏心 
距 , 两 端 为 相同 椭圆 轴承 支撑 , 则 根据 图 2 所 示 椭 圆 
轴承 建立 第 卡尔 坐标 系 ,X 为 水 平方 向 ,7 为 垂直 方 
向 ,Z 为 轴 向 方向 , 则 建立 转子 -轴承 系统 动力 学 方 


程 为 
闪避 Me cos(¥ | (22) 
My a Me(? sin(X - G 
式 中 :Ff、, 分别 为 轴承 对 转子 在 xy 方向 的 油膜 
力 , 对 上 式 进行 无 量 纲 化 ,可 得 无 量 纲 动力 学 方程 为 


Teosr 
网 (23) 
[| 


Rn 
圆 畏 承 在 非 承载 瓦 开 槽 及 改变 供 油 压力 时 刚性 转子 
的 动力 学 特性 与 稳定 性 。 


3. 计算 结果 与 分 析 


二 刚性 转子 质量 为 2M =200 kg, 转 子 质心 偏心 距 
cE8. 8 x 10 mm。 椭圆 轴承 参数 :转子 直径 D = 
100Dnm; 轴 承 半径 间 阶 C=0.2 mm; 上 下 轴瓦 包 角 均 
将 访 0°; 上 下 瓦 的 预 负荷 系数 为 0.2; 润 滑 油 动力 血 
度 胃 . = 0.029 8 Pa . s; 轴 瓦 长 径 比 A, =0.9, 以 轴承 
咏 区 坐标 (0,0) 为 初始 值 ,计算 结果 如 图 5 ~ 11 
jj 全 . 


= 
O 


无 量 纲 位 移 了 


-To 
=1:05 0:8. 0;6 04 0.2 010 0.2 :0:40:6- "0:8 1:0 
无 量 纲 位 移 X 
图 5 转速 为 7000 r/min, 普 通 轴承 与 
槽 轴承 的 转子 轴 心 轨迹 
Fig.5$ At the rotation speed of 7 000 r/min ,rotor orbits of 


ordinary bearing and bearing with different width grooves 
图 5 为 转速 为 7 000 r/min 时 , 覃 圆 轴 承 非 承载 
瓦 开 不 同 槽 宽度 但 不 通 压 力 油 时 对 转子 轴 心 轨迹 的 
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大 ,转子 的 轴 心 轨迹 也 随 之 增 大 ,说 明 开 槽 降低 了 系 
统 的 稳定 性 。 

图 6、 图 7 为 同一 开 槽 宽度 不 同 供 油 压 力 下 , 系 
统 转子 的 频谱 图 。 通 过 图 6 与 图 7 的 对 比 ,可 以 看 
出 同一 转速 7 000 r/min 下 , 开 槽 宽度 一 定时 , 槽 内 
共 油 压力 的 增加 可 以 有 效 降 低 转 子 振动 幅 值 ,抑制 
转子 半 频 振动 成 分 , 半 频 振动 幅 值 由 0.54 降 至 
0. 005 ,降低 了 99% 。 

通过 图 7 与 图 8 的 对 比 ,可 得 当 供 油 压 力 一 定 
时 , 开 槽 宽度 的 增加 能 使 转子 振幅 进一步 降低 , 半 频 
振动 的 幅 值 由 0.005 降 为 0, 转 子 振动 成 分 由 半 频 成 
分 转变 为 基 频 成 分 ,有 效 提高 了 系统 的 稳定 性 。 


0.6 


雁 


QO 05 0 L320 23 30 33 di0 5 530 
振动 阶 次 
图 6 转速 为 7000 r/min,L,=0.1L,p, =0.5 MPa 时 ， 
转子 振动 频谱 图 
Spectrum diagram of rotor vibration at 7 000 r/min, 
L,.=0.1L,p,=0.5 MPa 

0.035 
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0.025 


夯 
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0.015 


大约 拓 | 


0.010 
0.005 
0.000 
LO 0 1:0° 1:5. 2:0° .25112300 3540 45 5.0 
振动 阶 次 
图 7 转速 为 7000 r/min,L,=0.1L,p, =0.8 MPa 时 ， 
转子 振动 频谱 图 
Fig.7 Spectrum diagram of rotor vibration at 7 000 r/min, 
L,.=0.1L,p, =0.8 MPa 

图 9 ~ 图 11 分 别 为 不 同 转速 下 系统 转子 振动 频 

谱 图 。 图 9 为 不 开 槽 的 普通 轴承 的 情况 。 从 图 9 可 

以 看 到 转子 在 普通 椭圆 轴承 里 以 转速 5 300 r/min 
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运行 时 ,频谱 中 开始 出 现 半 频 振 动 成 分 。 
0.035 
0.030 
0.025 
] 吗 
媒 0.020 
以 
噶 0.015 
全 
0.010 
0.005 
0.000 
O00 05 有 0 下 5 人 5 
振动 阶 次 
图 8 转速 为 7000 r/min,L, =0.3L,p, =0.8 MPa 时 ， 
转子 振动 频谱 图 


Fig.8 Spectrum diagram of rotor vibration at 7 000 r/min, 
L.=0.3L,,p, =0.8 MPa 
本 由 图 10 和 图 11 可 知 ,改变 轴瓦 结构 和 供 油 压 
万 局 ,能 使 系统 出 现 半 频 振动 时 的 转速 明显 提升 , 提 
浪 包 系统 的 稳定 性 。 
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振动 阶 次 
转速 为 5 300 min ,转子 在 普通 轴承 里 振动 频谱 图 
9 The vibration spectrum of the rotor for 5 300 r/min 
0.040 


chigaXiv 
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振动 阶 次 
10 ”转速 为 8 400 r/min,L. =0.3L,p, =0.8 MPa 时 ， 
转子 振动 频谱 图 
Fig. 10 Spectrum diagram of rotor vibration at 8 400 r/min, 
L,=0.3L,p, =0.8 MPa 
图 10 和 图 11 为 开 槽 通 油 后 的 情况 , 开 权 宽 度 
为 0.37, 供 油 压 力 为 0.8 MPa。 从 图 10 可 知 ,频谱 


Es 
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图 中 显示 转子 转速 为 8 400 min 半 频 振动 刚 出 现 ， 
相 比 于 普通 椭圆 轴承 转速 提升 了 58.5% 。 从 图 11 
可 知 ,转子 转速 为 8 500 r/min 时 , 半 频 振动 振幅 进 
一 步 增 大 。 
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图 11 转速 为 8500 r/min,L, =0.3L,p, =0.8 MPa 时 ， 
转子 振动 频谱 图 
Fig.11 Spectrum diagram of rotor vibration at 8 500 r/min, 


L, =0. 3L,,p; =0.8 MPa 


本 研究 提出 根据 动态 5 油膜 的 边界 条 件 求解 油 
膜 负 压 与 轴 向 压 差 共同 作用 下 的 压力 分 布 的 方法 ， 
得 到 了 轴 向 压 差 非 零 时 椭圆 轴承 轴瓦 上 的 压力 分 布 
近似 解 , 解 决 了 工程 应 用 中 椭圆 轴承 非 承 载 区 周 问 
压力 槽 内 油膜 压力 非 线性 求解 问题 。 运 用 该 方法 快 
速 求解 得 到 了 椭圆 轴承 存在 轴 向 压 差 时 的 非 线性 雷 
庄 方 程 解 ,并 与 有 限 差 分 法 对 比 , 保 证 所 得 解 的 正 
确 性 。 

运用 所 提 方 法 ,研究 了 椭圆 轴承 非 承 载 区 周 向 
槽 内 压力 对 近似 解析 解 的 影响 ,得 到 适用 于 在 动 压 
椭圆 轴承 非 承载 瓦 开 覃 的 计算 模型 ,计算 后 得 到 如 
下 结论 。 

1) 椭圆 轴承 长 径 比 取 不 同 值 时 ,由 于 预 负 奏 系 
数 以 及 槽 内 供 油 压力 的 存在 ,轴承 对 转子 的 油膜 力 
在 一 定 范围 内 并 不 随 油 膜 力 的 增 大 而 增 大 。 

2) 在 非 承 载 上 瓦 只 开 圆 周 槽 不 通 压力 油 会 使 
系统 的 稳定 性 降低 。 

3) 非 承载 上 瓦 开 槽 并 供 油 时 ， 开 槽 党 度 一 定 
时 ， 增 大 覃 内 供 油 压 力 能 降低 转子 振动 幅 值 ， 抑 制 
转子 振动 的 半 频 成 分 ; 在 供 油 压力 一 定时 ， 增 大 开 
槽 宽度 ， 能 进一步 抑制 转子 振动 幅 值 ， 提 高 系统 稳 
定性 。 
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